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Abstract 

We consider the g-state Potts model on families of self-dual strip graphs of the square 
lattice of width Ly and arbitrarily great length L^., with periodic longitudinal boundary 
conditions. The general partition function Z and the T = antiferromagnetic special case P 
(chromatic polynomial) have the respective forms Y^^=i^'^ CF,Ly,j{^F,Ly,j)^'' , with F = Z,P. 
For arbitrary Ly, we determine (i) the general coefficient CF,Ly,j in terms of Chebyshev 
polynomials, (ii) the number nF{Ly,d) of terms with each type of coefficient, and (iii) the 
total number of terms Np^Ly^x- We point out interesting connections between the nz{Ly,d) 
and Temperley-Lieb algebras, and between the Np^^^^x and enumerations of directed lattice 
animals. Exact calculations of P are presented for 2 < Ly < 4. In the limit of infinite length, 
we calculate the ground state degeneracy per site (exponent of the ground state entropy), 
W{q). Generalizing q from Z_|_ to C, we determine the continuous locus B in the complex q 
plane where W{q) is singular. We find the interesting result that for all Ly values considered, 
the maximal point at which B crosses the real q axis, denoted qc is the same, and is equal 
to the value for the infinite square lattice, qc = 3. This is the first family of strip graphs of 
which we are aware that exhibits this type of universality of qc- 
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1 Introduction 



The g-state Potts antiferromagnet (AF) [|l], |^ exhibits nonzero ground state entropy, 5*0 > 
(without frustration) for sufficiently large g on a given lattice A or, more generally, on a 
graph G = (V, E) defined by its set of vertices V and edges joining these vertices E. This 
is equivalent to a ground state degeneracy per site W > 1, since Sq = kslnW. There is 
a close connection with graph theory here, since the zero-temperature partition function of 
the above-mentioned g-state Potts antiferromagnet on a graph G satisfies 

Z{G,q,T = 0)pAF = P{G,q) (1.1) 

where P{G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
(for reviews, see @-0). The minimum number of colors necessary for such a coloring of G 
is called the chromatic number, x(G'). Thus 

W{{G},q)= hm P(G,g)i/" (1.2) 

71— >00 

where n = \V\ is the number of vertices of G and {G} = lim„^oo G. At certain special points 
qs (typically g^ = 0, 1, .., x{G)), one has the noncommutativity of limits 

lim lim P(G,g)^/" ^ lim lim P(G,g)i/" (1.3) 

q—*qs n—*oo n— >oo q—*qs 

and hence it is necessary to specify the order of the limits in the definition of PF({G}, g^) [0. 
Denoting Wgn and Wnq as the functions defined by the different order of limits on the left 
and right-hand sides of ( p..3| ), we take W = Wgn here; this has the advantage of removing 
certain isolated discontinuities that are present in Wnq- 

Using the expression for P{G,q), one can generalize g from Z+ to C The zeros of 
P(G, g) in the complex g plane are called chromatic zeros; a subset of these may form an 



accumulation set in the n — oo limit, denoted B [^, which is the continuous locus of points 
where W{{G},q) is nonanalytic. (For some families of graphs B may be null, and W may 
also be nonanalytic at certain discrete points; this is not relevant for the present paper.) The 
maximal region in the complex g plane to which one can analytically continue the function 
VrdG}, g) from physical values where there is nonzero ground state entropy is denoted Pi. 
The maximal value of g where B intersects the (positive) real axis is labelled qc{{G}). This 
point is important since it separates the interval of physical g > qc{{G}) on the positive real 
g axis where the Potts model exhibits nonzero ground state entropy (without frustration) 
from the interval < g < qc{{G}) in which W has different analytic form(s). 
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In this paper we shall present a number of exact results on the Potts model partition 
function and the special case comprised by the zero-temperature Potts antiferromagnet for 
families of self-dual strip graphs of the square lattice with (i) a fixed transverse width Ly, 
(ii) arbitrarily great length L^, (iii) periodic longitudinal boundary conditions, and (iv) such 
that each vertex on one side of the strip, which we take to be the upper side (with the strip 
oriented so that the longitudinal, x direction is horizontal) are joined by edges to a single 
external vertex. A strip graph of this type will be denoted generically a,s Gd (where the 
subscript D refers to the self-duality) and, when its size is indicated, as Gn^Ly x L^)- We 
shall present a number of structural formulas that hold for arbitrary Ly. After recalling 
our earlier work for the case Ly = 1 0, 0, we shall present exact results for the chromatic 
polynomial for Ly = 2 through Ly = 4 and shall study the zeros of these polynomials and 
determine the continuous accumulation set B of these zeros in the limit L^ oo. In Fig. 
ID we show an illustrative example of this family of strip graphs, for the case L^ = 4 and 
L,, = 3. 
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Figure 1: Example of a strip graph GoiLy x L^) for the 



In general, the graph GoiLy x L^) has 

n=\V\=L,Ly + l = f (1.4) 

vertices, equal to the number of faces, /, and 

e = \E\= 2L^Ly (1.5) 

edges. We recall that the dual G* of a planar graph G = (V, E) is defined as the graph 
obtained by replacing each vertex (face) of G by a face (vertex) of G* and connecting the 
vertices of the resultant G* by edges (bonds). The graph G is self-dual if and only ii G = G*. 
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It is easily checked that the family GoiLy x L^) is self-dual. We recall that the degree A of a 
vertex is defined as the number of vertices to which it is adjacent (i.e., the number of nearest 
neighbors). In the graph GoiLy x Lx) there are Lx{Ly — 1) vertices with degree A = 4, 
Lj. vertices with A = 3, and finally the one external vertex with A = L^. The chromatic 
number is 

The (infinite) square lattice (interpreted as the infinite limit of a graph) is self-dual, and 
this property has been useful in studies of the Ising, and general g-state Potts model on this 
lattice. An early example of this was the calculation by Kramers and Wannier of the critical 
temperature of the (zero-field) Ising model on the square lattice [§] before Onsager derived 
a closed-form expression for the free energy ^ . The Potts model partition function and the 



equivalent Tutte polynomial fl^, |Tl| of a planar graph G satisfy certain symmetry relations 



when one changes G to G* p|, |T^. Let us define v = — 1, where K = J/{kBT), with T 



being the temperature and J the spin-spin exchange constant for the Potts model. Then 0] 

Z{G,q,v) = v<''\-<'''>Z{G*,qA) . (1.7) 

V 

li G = G*, this reduces to an identity relating the g-state Potts model on G to the same 
model with v q/v. In previous studies of complex-temperature (Fisher) zeros of the 
Potts model, it has proved useful to employ families of sections of the square lattice that 
are self-dual so that one can take advantage of the symmetry property (|1.7| ) |ll3l-|l^. The 
resulting patterns of zeros can be compared with those obtained with other types of boundary 
conditions |l^, [l^, |20|. It may be recalled that in finite-lattice studies one often chooses 
periodic boundary conditions to minimize boundary effects, but these do not preserve the self- 
duality property of the infinite square lattice. There are several types of self-dual boundary 
conditions for strips of the square lattice Jl3|, |l^, [l^. (The boundary conditions involved 
in the family Go^Ly x L^) were denoted DBC2 in [|^.) In this paper, instead of studying 
complex-temperature zeros, we shall focus on zeros of P{G, q) in the complex q plane. As we 
shall show, these exhibit some very interesting features. In particular, we shall show that for 
all of the strip widths considered, we find a universal value of qc, namely qc = 3, the value 



for the infinite square lattice 121 



2 Some General Properties 

For a recursive family of graphs, such as the strip graphs considered in this paper, comprised 
of repetitions of a basic subgraph, the g-state Potts model partition function for arbitrary 
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q (not necessarily in Z+) and v has the form |jT2 



Z{G,q,v)= Cz,G,A^z,G,f^ (2.1) 

3=1 

where the coefficients cz,Gj and the terms Xz,G,j depend on the lattice type, the boundary 
conditions, and the width, but not the length. Since we are only dealing in this paper with 
a particular type of strip graph, Gr,, we shall henceforth usually suppress the Gd in the 
notation, but shall make explicit the dependence on Ly, writing Nz^Ly,x for the total number 
of terms, Xz,Ly,j for a given term, cz,Ly,j for the corresponding coefficient, and similarly for 
the chromatic polynomial. Since the chromatic polynomial is a special case of the Potts 
model partition function, it also has this form |2r 



^P,Ly,X 

P{Gn{LyXL,),q)= ^ cp,^„,(Ap,i^,,)^^ . (2.2) 

3=1 

Following our earlier nomenclature ||^, we denote a A as leading (= dominant) if it has a 
magnitude greater than or equal to the magnitude of other A's. In the limit n ^ oo the 
leading A in Z determines the free energy per site / = lim„^oo Z^^"" and similarly the leading 
A in P determines the function W defined in ( |1.2| ). The continuous locus B where f or W is 
nonanalytic thus occurs where there is a switching of dominant A's in Z and P, respectively, 
and is the solution of the equation of degeneracy in magnitude of these dominant A's. 

In general one can regard the A's as the (nonzero) eigenvalues of a certain coloring matrix 
2^ . In our earlier calculations, we have found examples of zero eigenvalues , but 



for the present family we only obtain nonzero eigenvalues. We note that while the coefficients 
cp,Lyj and cz,Lyj can be obtained as multiplicities of the distinct eigenvalues Xp,Ly,j and 
Xz,Ly,j for sufficiently large integer q, they may be zero or negative for some positive integer 
values of q. More generally, while they play the role of eigenvalue multiplicities for sufficiently 
large integer q, the domain of their definition may be generalized to g G M+ or, indeed, to 
qeC. 

The dimension of the space of coloring configurations, Af, is equal to the sum of the 
multiplicities of each distinct eigenvalue, i.e., the sum of the dimensions of the invariant 
subspaces corresponding to each of these distinct eigenvalues. For the chromatic polynomial, 
this is A/", which is equal to the sum 



Cp,Ly= Yl ^P^^y^^ (2.3) 
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while for the full Potts model partition function, we shall denote it as 

Cz,Ly = 5Z ^^'^y^ ■ 



(2.4) 



Some previous literature on chromatic polynomials and Potts model partition functions on 
lattice strips is in p7[|-[p8[ . 



3 Structural Theorems 



3.1 Coefficients 



We find that the coefficients that enter into ( ^^ for P{Gd,(i) and (|2.1|) for Z{GD,q,v) are 
polynomials in q that consist of a special restricted set with the property that the coefficient 
of maximal degree ci > 1 in g is 



id) 



u. 



2d 



V 2 



2d 



\ 2 



(^ 



d-l 
j=0 



(3.1.1) 



where Tn{x) and Un{x) are the Chebyshev polynomials of the first and second kinds, defined 
by 



TJx) 



n I n — J 



j=0 



n-J \ 3 



(2x) 



n-2j 



and 



i=o 



(3.1.2) 



(3.1.3) 



where in eqs. ( |3.1.2| ) and ( |3.1.2| ) [|] in the upper limit on the summand means the integral 
part of |. The first few of these coefficients are 



^(2) = g(g _ 2) 



^(3) = l)(g_3) 



(3.1.4) 
(3.1.5) 

(3.1.6) 
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=g(g_2)(g2_4g + 2) (3.1.7) 

t^^^) = qi^q^ - ^ l)i^q^ - ^ ^) (3.1.8) 

^(6) = g(g - l)(g - 2)(g - 3)(g2 - 4g + 1) . (3.1.9) 
We have established the following factorization: 

d 

^^'^ = U^q - sd,k) (3.1.10) 

k=l 

where 



Sd,k = 2 + 2 cos 



= 4cos2(^^) for A; = 1,2,.. d (3.1.11) 

In |]2^, we showed that the coefficients cl^j that entered into Z{Gs,q,v) and hence also 
P{Gs, q) for cychc strips G^of the square or triangular lattices (and Mobius strips of the 
square lattice) were of the form 



d 



'2d -3 



d 
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l[{q-qd,k) (3.1.12) 



k=l 

where 

gd,fc = 2 + 2cos(^^^) =4cos'(^^^) for A; = 1, 2, ...d . (3.1.13) 
The K^'^^ can be expressed as differences of these c^'^^ coefficients: 

^{d) ^ ^(d) _ ^{d-i) for ci= 1,2,... (3.1.14) 

The K^"^^ coefficients have the following general properties. 

• Given the expression of k^'^\ ( |3.1.1| ) in terms of Chebyshev polynomials, and the re- 
cursion relation for these Chebyshev polynomials (for n > 1) 

Chn+i{x) = 2xChn{x) — Chn-i{x) for Chn{x) = Tn{x) or Un{x) (3.1.15) 

we obtain the recursion relation for the k^'^^ for d > 2, 

^(d+i) = _ 2)k^'^^ - . (3.1.16) 
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• K^'^'^ is a polynomial of degree d in g whose coefficients alternate in sign. 

• k!^'^^ has the factor q. 

• If and only if d is even and d >2, then k*^'^-' has the factor (g — 2). 

• If and only if = mod 3 and d > 3, then n^'^^ has the factor (g — 3). 

• The coefficient of the term q'^ in k^'^'' is 1. 

• The coefficient of the term q in k^'^^ is (— 

• The d zeros of k^'^'' occur at the (real) values SdM, k = 1, ...,d, where s^.fc was given in 
eq. ( 3.1.11 ). Clearly, these lie in the interval < g < 4. In passing, we note that the 
k = 2 special case 2 is the Tutte-Beraha number Bd- 

If one lets 

g = 2 + 2cos^ = 4cos2(^-j (3.1.17) 

one sees that the argument of the Chebyshev polynomials in ( ^.1.1| ) is given by ^ = cos(0/2). 
Using the identities (with uj = 9/2 here) 

T„(cosct;) = cos{nuj) (3.1.18) 

and 

TT ( \ sin[(ra + l)w] 

t/„(cosc<j) = — (3.1.19) 

sinu; 

we have an alternate formula for k'^'^-', 

K^'^^ = 2 cot sin(2dcj) (3.1.20) 



3.2 Determination of np{Ly,d) and Np^i^^x 

Let us define np{Ly,d) as the number of terms \p,Ly,j in PiGoiLy x L^),g) that have as 
their coefficients cp^^^j = k^'^\ For G£>{Ly x L^) strips, these coefficients are nonzero for 
1 < d < Ly + 1. (In our earlier work |2^ for different strip graphs, the coefficients were 
defined to start with S^'^ = 1, so that the range went from < (i < Ly] in both cases, the 
resultant chromatic polynomial has terms in g ranging from g" to g.) The total number, 
Np,Ly,\, of different terms \p,Ly,j in ( p.2| ) is then given by 

Ly + l 

Np^Ly,x=J2npiLy,d). (3.2.1) 

d=l 
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For the sum, ( p.3|) , of the coefficients in (|2.2| ) we have 



Ly + 1 



(3.2.2) 



d=l 



Using the same methods as in ||25|, we have determined the numbers np{Ly, d) of Xp^Ly/s 
that have each type of coefficient, k^'^^ From coloring matrix arguments, using the fact that 
a transverse shce along the strip is a path (tree) graph T^^+i, for which the chromatic 
polynomial is P(Tiy+i, q) = q{q — 1)^", it follows that 



Ly+l 



J2 MLy, d)K^^^ = P{TLy+i, q) = q{q - 1)' 



(3.2.3) 



Differentiating eq. ( p.2.3| ) Ly times we obtain Ly + 1 linear equations in the Ly + 1 unknowns 
np{Ly, d) for 1 < d < Ly + 1. Solving these equations, we find that 



np{Ly, d) = for d > Ly + 1 
np{Ly,Ly + 1) = 1 



and 



(3.2.4) 
(3.2.5) 

(3.2.6) 



np(l,l) = l 

with all other numbers np{Ly,d) being determined by the recursion relation 

np{Ly + l,d) = np{Ly,d -1) +np{Ly,d) + np{Ly,d+l) for 1 < d < Ly + 1 . (3.2.7) 
In particular, we find 



np{Ly,Ly) - Ly 



and 



np{Ly, 1) = MLy 
where M„ is the Motzkin number, given by 



(3.2.8) 
(3.2.9) 

(3.2.10) 
(3.2.11) 

is the Catalan number. (The Catalan and Motzkin numbers occur in many combinatoric 
applications ||59||-||6l||.) 



i=o ^-^ ^ 



and 



1 (In 
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Table 1: Table of numbers np(Ly,d) and their sums, Np,L^^\ for GoiLy x Lx). Blank entries are zero. 



Ly I d ^ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 




1 


1 


1 




















2 


2 


2 


2 


1 


















5 


3 


4 


5 


3 


1 
















13 


4 


9 


12 


9 


4 


1 














35 


5 


21 


30 


25 


14 


5 


1 












96 


6 


51 


76 


69 


44 


20 


6 


1 










267 


7 


127 


196 


189 


133 


70 


27 


7 


1 








750 


8 


323 


512 


518 


392 


230 


104 


35 


8 


1 






2123 


9 


835 


1353 


1422 


1140 


726 


369 


147 


44 


9 


1 




6046 


10 


2188 


3610 


3915 


3288 


2235 


1242 


560 


200 


54 


10 


1 


17303 



We find that the np{Ly,d) are closely related to certain numbers enumerating random 
walks. To explain this, consider a random walk on the nonnegative integers such that in 
each step the walker moves by +1, —1, or units. Denote m{n, k) as the number of walks of 
length n steps starting at and ending at k. The case k = describes the number of walks 
defined above that return to the origin after n steps. This is given by the Motzkin number; 
m{n, 0) = Mn- Define the sum 

n 

5„ = ^m(n,A;). (3.2.12) 

A!=0 

We find that the numbers np{Ly,d) are precisely equal to these enumerations of random 
walks: 

np{Ly,d) =m{Ly,d-l) ioi 1 < d < Ly + 1 (3.2.13) 

The values of the np{Ly, d) for Ly = 1 to Ly = 10 are listed in Table |l|. (Note that Sg should 
read 6046 in the related Table 3 of Ref. 

From the identity (|3.2.13| ), it follows that the total number of terms Np^i^^x is 

Np,L„x = Sl, . (3.2.14) 

These numbers are listed in Table ^ 

We have determined a generating function whose expansion yields the numbers Np^i^y. 



2 



{^) - l] - X = Y: NP,L,,^-'^^' . (3.2.15) 
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Combining this with our results from [^], we observe the interesting fact that the total 
numbers of terms Np^^^^x ioT the chromatic polynomial of a G^iLy x L^) strip graph (for 
any L^) is equal to 1/2 the corresponding total number of terms for the cyclic or Mobius 
strips of the square or triangular lattice with the next larger width, Ly + 1, i.e., 



Np,Ly,\ = -Np^A,cyc/Mb,Ly+i A = sq,tri . 



Further, we observe the intriguing relation 



DA,sq,Ly+l 



(3.2.16) 



(3.2.17) 



where Ni:)A,sq,n denotes the number of directed n-site lattice animals on the square lattice. 
We prove this by noting that the generating function (|3.2.15| ) that we have found for Np^L^^x 
is the same as the generating function for directed n-site lattice animals on the square lattice 
given in [|6^ . 

From ( |3.2.15| ), it follows that for large width L^, Np^i^^x grows exponentially fast, with 
the leading asymptotic behavior: 



AT, 



p,L„A ~ Ly^'^ 3^« as Lj, ^ oo . 



(3.2.18) 



This is the same asymptotic behavior that we found in for the total number of terms en- 
tering in the chromatic polynomials for the cyclic or Mobius strips of the square or triangular 
lattice of width L,, 



3.3 Determination of nz{Ly,d) and Nz^Ly,x 

From the same type of coloring matrix arguments used in we have 



Ly + l 



(3.3.1) 



d=l 



Differentiating eq. ( p. 3.11) Ly times we obtain Ly + l linear equations in the Ly + l unknowns 
nz{Ly, d) for 1 < d < Ly + 1. Solving these equations, we obtain the general results 



and 



uziLy, d) = for d > Ly + 1 
nz{Ly,Ly + 1) = 1 

np(l,l) = 2 



(3.3.2) 
(3.3.3) 

(3.3.4) 
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with all other numbers nz{Ly, d) being determined by the recursion relation 
nz{Ly + 1, d) = nz{Ly, d — 1) + 2nz{Ly, d) + nz{Ly, d+ 1) for 1 < d < Ly + 1 . (3.3.5) 



We solve this recursion relation with the conditions (|3.3.2| )- (|3.3.4| ) in closed form and obtain 



Summing these numbers, we find, for the total number of terms 

N,.., ^ {fji) . (3.3.7) 

We observe that the total numbers of terms Nz,Ly,\ for the Potts model partition function 
on the Go^Ly x L^) strip graph (for any L^) is equal to 1/2 the corresponding total number 
of terms for the cyclic or Mobius strips of the square or triangular lattice with the next larger 
width, Ly + 1, i.e., 

Nz,Ly,\ = -Nz,A,cyc/Mb,Ly+i A = sq,tri . (3.3.8) 
There is also an interesting connection with directed lattice animals: 

Nz,Ly,X = NDA,tri,Ly+l (3.3.9) 

where Nj:)A,tri,n denotes the number of directed ra-site lattice animals on the triangular lattice. 
As — > oo, Nz,Ly,\ has the leading asymptotic behavior 

Nz,Ly,X ~ 7r-l/2^,;l/2 ^Ly Ly ^ oo . (3.3.10) 

The values of nz{Ly, d) and Nz,Ly,d are given for 1 < L^, < 10 in Table ||. 

In what follows, we shall need to refer to our calculation of the numbers nz{Ly, d) of A's 
with associated coefficient c^^^ for cyclic strips Gs of the square or triangular lattice of width 



Ly (and arbitrarily L^ in ||2^. We found 



for < d < Ly and zero otherwise. These values will be needed below and are listed in Table 
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Table 2: Table of numbers nz(Ly,d) and their sums, Nz,Ly.\ for GniLy x Lx). Blank entries are zero. 



Ly I d ^ 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


Np,Ly,X 


1 


2 


1 




















3 


2 


5 


4 


1 


















10 


3 


14 


14 


6 


1 
















35 


4 


42 


48 


27 


8 


1 














126 


5 


132 


165 


110 


44 


10 


1 












462 


6 


429 


572 


429 


208 


65 


12 


1 










1716 


7 


1430 


2002 


1638 


910 


350 


90 


14 


1 








6435 


8 


4862 


7072 


6188 


3808 


1700 


544 


119 


16 


1 






24310 


9 


16796 


25194 


23256 


15504 


7752 


2907 


798 


152 


18 


1 




92378 


10 


58786 


90440 


87210 


62016 


33915 


14364 


4655 


1120 


189 


20 


1 


352716 



Table 3: Table of numbers nz{Ly, d) and their sums, Nz,g,x for cyclic strips Gs{Ly x L^.) of the square or 
triangular lattice of width Ly and arbitrary length . 



Ly [ d ^ 





1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


Nz,Ly,X 


1 


1 


1 




















2 


2 


2 


3 


1 


















6 


3 


5 


9 


5 


1 
















20 


4 


14 


28 


20 


7 


1 














70 


5 


42 


90 


75 


35 


9 


1 












252 


6 


132 


297 


275 


154 


54 


11 


1 










924 


7 


429 


1001 


1001 


637 


273 


77 


13 


1 








3432 


8 


1430 


3432 


3640 


2548 


1260 


440 


104 


15 


1 






12870 


9 


4862 


11934 


13260 


9996 


5508 


2244 


663 


135 


17 


1 




48620 


10 


16796 


41990 


48450 


38760 


23256 


10659 


3705 


950 


170 


19 


1 


184756 
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Table 4: Generic Bratteli diagram relevant for TLn{q). 



n [ m 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


11 


1 




1 




















2 


1 




1 


















3 




2 




1 
















4 


2 




3 




1 














5 




5 




4 




1 












6 


5 




9 




5 




1 










7 




14 




14 




6 




1 








8 


14 




28 




20 




7 




1 






9 




42 




48 




27 




8 




1 




10 


42 




90 




75 




35 




9 




1 



4 Connection with Bratteli Diagrams for TLn{q) 



There are a number of interesting algebraic properties of structural elements and properties 
of the chromatic polynomials and Potts model partition functions equivalently Tutte polyno- 
mials. The interpretation of the coefficients cp^Csj chromatic polynomials as dimensions 
of invariant subspaces of coloring matrices was discussed in [E^, |23] and its generalization to 



the full Potts model partition function was discussed in p5|. Thus, Cp^Gs,Ly and Cz,Gs,Ly are 
sums of these dimensions, a fact that we have used in and eqs. (|3.2.3| ) and (|3.3.1|) here. 



In 1 25] we pointed out another connection, namely the fact that the numbers nz{Gs, Ly, d) 



for < d < Ly for cyclic strips Gg of the square or triangular lattice of width Ly were 
related to Temperley-Lieb algebras. We have found a very interesting generalization of this 
for the nziGo, Ly,d) calculated here. Let us denote TLn{q) as the Temperley-Lieb algebra 
of operators Ui, 



n 



, n satisfying |6^, |6^ 



U 



[f/„f/,]=0 if^-jl^l 



(4.1) 
(4.2) 



where \A, B\= AB - BA, and 



[/,f/,±if/, = U, (4.3) 
Associated with the structural decomposition of this algebra TLn{q) there is what is known 



Bratteli diagram ||6^, ^ , given in Table |. 
We observe the following interesting properties: 
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• For even n = 2Ly, the entries in this row of the Bratteh diagram are equal to the 
nz{Ly, d) for the cychc strip of the square lattice of width Ly and arbitrary length L^. 
for m = 2d + 1 with d in the range < d < Ly, i.e., 1 < m < 2Ly + 1 = n + 1. 

• For odd n = 2Ly + 1, the entries in this row of the Bratteli diagram are equal to the 
nz{Ly, d) for the strip Gd of width Ly and arbitrary length L^ for m = 2d with d in 
the range I < d < Ly + 1, i.e., 2 < m < 2{Ly + 1) = n + 1. 

As this makes clear, the numbers nz{Ly,d) for the cyclic strip and the Gd strip can be 
expressed via a single formula. Define 

Q».=r":'Vf!"l) (4.4) 



m J \m — 2^ 
so that the Cn,m satisfy the recursion relation 

Cn,m = Cn-l,m-l + C*n-l,m+l • (4-5) 



For Gcyc{Ly X Lx), let n = 2Ly and m = Ly — d with < d < Ly. For Go^Ly x L^), let 
n = 2Ly + 1 and m = Ly + 1 — d with 1 < d < Ly + 1. Then with these substitutions, C, 
reduces to the respective expressions ( |3.3.11|) and ( p.3.6| ). 



Indeed, the fact that our results in related nz{Ly,d) for the cyclic strips of the 
square and triangular lattices to the even-ra rows of the Bratteli diagram, i.e., the rows for 
TL2Ly{(l), raised a question: is there a family of graphs with the property that the coefficients 



in the generic formula (2J.) are restricted to a simple set such that the number of A's with 
coefficients of a particular type fill out the odd-n rows of the Bratteli diagram. We had 
searched for such a family of graphs and had finally found it in the form of the GoiLy x Lx) 
graphs, which have provided an affirmative answer to this question: the entries in the odd 
n = 2Ly + 1 rows are given by the nz{Ly, d), as described above. Thus, our current results 
are of interest both for the connection with Bratteli diagrams and for the fact that, together 
with our previous findings in |^5[, they saturate this connection, i.e., the entries in all of the 
Bratteli diagram are now accounted for in terms of the 1-1 correspondence with nz{Ly,d) 
numbers for these cyclic and Gd strip graphs. 

We observe that the Bratteli diagram also has a combinatoric interpretation. Consider 
random walks on the non-negative integers that begin at the origin and have the property 
that the walker moves by -|-1 or —1 spaces at each step. The (n, m) entry in the Bratteli 
diagram is the number of random walks of this type, consisting of n steps, that end at m — 1. 
Note that the number of random walks of this type that consist of 2n steps (a misprint in 



had this as n steps) is the Catalan number C„. The recursion relation ( [4.5|) follows 
immediately from this. 
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5 Width Ly = 1 

It is of interest to calculate Potts model partition functions and chromatic polynomials for 
these families of self-dual graphs. In this section we shall discuss the chromatic polynomial 
for the lowest case, Ly = 1, and in the next sections we shall present calculations of the 
chromatic polynomials for 2 < Ly < 4 and shall analyze the zeros of these polynomials in 
the complex q plane and their accumulation sets as — cx). In a separate publication we 
shall present more lengthy calculations of the full Potts model partition functions. 

We recall that the Ly = 1 family is comprised of the wheel graphs (Wh) with L^ vertices 
on the rim and a central vertex, with edges forming "spokes" connecting this central graph 
to the vertices on the rim. The n-vertex wheel graph can be constructed as the "join" 
(Wh)n+i = Ki + Cn, where Kn and C„ are the complete graph and the circuit graph on 
n vertices, respectively. (The join of two graphs G and H is defined as the graph obtained 
by connecting each vertex of G to every vertex of H by edges. The complete graph Kn is 
defined as the graph with n vertices such that each vertex is adjacent to every other vertex.) 
The chromatic polynomial has the form of (|2.2|) with Np^ =i a = 2 and the terms 



Ai,i = -1 (5.1) 

Ai,2 = g-2 (5.2) 

and coefficients 

ci,i = K^^^ (5.3) 

ci,2 = . (5.4) 

Here and below for brevity we write XLyj for Xp^Lyj and CLyj for cpxyj- Thus, 

P(l X L,,g) = {Wh)L^^, = q[iq- 2)(-l)^^ + (g - 2)^^] . (5.5) 

The continuous locus B formed as the continuous accumulation set of chromatic zeros in the 
limit La; oo hmit, is the unit circle centered at g = 2 ||^, |^ 

B: \q-2\ = l ioi Ly = l . (5.6) 

Evidently, this intersects the real axis at g = 1 and q = 3. The locus B separates the q 
plane into two regions, Ri and R2, which are, respectively, the exterior and the interior of 



the circle ( |5.6| ). In region Ri 

W = q-2 for g G i?i . (5.7) 
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In any other region than Ri, only the magnitude \W\ can be determined unambiguously p, 
and we have 

\W\ = 1 for g G i?2 • (5.8) 

Except for the single discrete zero at g = 2, the chromatic zeros lie exactly on the circle 
|g — 2| = 1 and occur with a density that is a constant as a function of angular position on 
this circle 0. 



6 Width Ly = 2 

Our general formula ( p.2.15| ) yields the result that there are A^p,Lj,=2,a = 5 terms for this 
width Ly = 2. The chromatic polynomial has the form (|2.2| ) with the following terms, listed 
in order of descending degree of their associated coefficients. 



A 



2,(2,3) 



A2,i = 1 (6.1) 
-(5-2g±y5) (6.2) 



A2,(4,5) - 2 L 



2 

- 5g + 7 ± (g^ - 6g^ + 15g^ - 22g + Uf^^ (6.3) 
and the corresponding coefficients 

C2,i = f^^^^ (6.4) 

C2j = k(2) forj = 2,3 (6.5) 

C2,j = K^^) for j = 4, 5 . (6.6) 

The locus B is shown in Fig. H, together with chromatic zeros for a long finite strip. 
Evidently, the chromatic zeros (except the discrete zero at g = 2), which eventually merge 
to form this locus in the oo limit, generally lie close to the curves on B. This locus 

divides the g plane into three regions: (i) region Ri including the intervals g > 3 and g < 1 
on the real axis and extending outward to infinitely large |g|; (ii) the region R2 including 
the real interval 5/2 < g < 3, and (iii) the innermost region R3 including the real interval 
1 < g < 5/2. There are triple points at which all three regions are contiguous at g = gt, g^, 
where qt — 5/2 + 1.27i. Note that the part of the boundary B separating regions R2 and R3 
is the line segment Re{q) = 5/2, —Im{qt) < Im{q) < Im{qt). In region Ri, A2,4 is dominant, 
so 

W = (A2,4)^/^ for qeRi (6.7) 
= |A2,3r/^ forgGi?2 (6.8) 
\W\ = \X2,2\^^^ for qeRs. (6.9) 
As is evident in Fig. ||, the locus B crosses the real g axis at g = 1, g = 5/2, and g = 3. 
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I. 



Re(q) 



Figure 2: Singular locus B for the ^ og limit of Gd(2 x L^. For comparison, chromatic zeros are 
shown for = 30 (i.e., n = 61). 
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7 Width Ly = ?> 

Here our general formula ( p.2.15|) yields Np^^^^^ x = 13. We have calculated the chromatic 
polynomial and obtain the following results, where again the A's are listed in order of de- 
scending degree of their coefficients. The first term is 

A3,i = -1 . (7.1) 

Next, the X^j for 2 < j < 4 are the roots of the equation 

- (3g - 8)^2 + (3g2 - 16q + 19)^ - (g^ - 8q^ + 19q - 13) = . (7.2) 

The Asj- for 5 < j < 9 are the roots of the equation 

+ (3g2 - 16g + 23)^^ + (3g^ - 34g3 + U2q^ - 261g + 177)^^ 

+ (g^ - 20q^ + 153g^ - 591g^ + 1234g2 - I332g + 581)^^ 

- (2g^ - 35g^ + 257g^ - 1028g^ + 2423g3 - 3371g^ + 2566g - 823)^ 

+ (g-2)2(g6_i4g5^79g4_ 230g3 + 366g2- 304g + 103) = . (7.3) 
Finally, the X^j for 10 < j < 13 are the roots of the equation 

- {q^ - 8g2 + 24g - 26)^^ - (2g^ - 23g^ + lllg^ - 279g^ + 362g - 191)^^ 

- {q- 3){q^ - Uq^ + 81g^ - 251g^ + 441^^ - 415g + 161)^ 

+ (g-2)3(g^-llg^ + 48g^-104g2 + 112g-47) = . (7.4) 
The corresponding coefficients are 

C3,i = (7.5) 

= k(^) for 2 < j < 4 (7.6) 

C3J = for 5 < j < 9 (7.7) 

C3,j = K^^^ for 10 < J < 13 . (7.8) 

The locus B defined in the Lx —>■ oo limit is shown in Fig. ^ together with chromatic 
zeros for a long finite strip. This locus divides the q plane several regions, including (i) Ri, 
which contains the real intervals q < 1 and g > 3 and extends outwards to infinite |g|; (ii) 
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Re(q) 



Figure 3: Singular locus B for the ^ oo limit of Gd(3 y. L^)- For comparison, chromatic zeros are 
shown for = 30 (i.e., n = 91). 
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i?2, which contains the real interval from q ~ 2.637 to g = 3; (iii) -R3, which contains the real 
interval 1 < g < 2.637; and (iv) the complex-conjugate pair of small regions i?4, centered 
approximately at g ^ 0.44 ± 1.672. The density of chromatic zeros varies in different parts 
of B] this density is low on the portion of B that passes through g = 1 and on the inner 
sides of the small regions R4^,R\. In the dominant term is a root of the quartic equation 
( [7.4|) , which we denote As^j^i, so that 

W = (A3,m)'/' for g e i?i . (7.9) 

In regions R2 and i?3, the dominant terms are two different roots of the fifth-degree equation 
([7. 3D, and a root of this equation is also dominant in the complex-conjugate pairs of regions 
i?4, R\. There are three complex-conjugate pairs of triple points on B. Our experience with 
the Lx ^ 00 limits of other families of strip graphs showed that there can often be additional 
extremely small regions [O; we have not made an exhaustive search for these. 



8 Ly = A 



Here our general formula (|3.2.15|) gives Np^4^x = 35. We have 

A4,i = 1 (8.1) 

with coefficient 

C4,i = K^^'> (8.2) 

and 

A4,2 = 2 - g . (8.3) 
The A4J for 3 < j < 5 are roots of the equation 

e + 3(g - 3)e - 3(g - 2)(g - 4)^ + (g=^ - Og^ + 24g - 17) = . (8.4) 

These X4J, 2 < j < 5 have the coefficient 

C4j = k(^) for 2 < j < 5 . (8.5) 

The other 30 terms consist of roots of two equations of degree 9 and an equation of degree 
12. Since these are rather lengthy, they are relegated to the appendix of the copy of this 
paper in the cond-mat archive. As before for the other Gd strips studied in this paper, these 
equations are in 1-1 correspondence with the coefficients listed in Table 0, so that all of the 
roots of a given equation have the same coefficient, as specified in this table. 
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The locus B defined in the Lx — * oo limit is shown in Fig. ^ together with chromatic 
zeros for a long finite strip. This locus B divides the q plane into a number of regions. These 
include (i) Ri, which contains the real intervals g < 1 and g > 3 and extends outwards to 
infinite |g| as before; (ii) R2, which contains the real interval from 2.746 ^ g < 3; (iii) a 
narrow region i?3, which contains the real interval (l/2)(3 + v^) ^ g ^ 2.746, and (iv) R^, 
which contains the real interval 1 < g ^ (l/2)(3 + -\/5); (v) the complex-conjugate pair of 
small regions R5, R^, centered approximately at g ~ ±1.4i; and (vi) a complex-conjugate 
pair of tiny sliver regions Rq, R^ located around g ~ 5/2 ± 1.55i. We note in passing that 
(l/2)(3 + = B5 = 2.61803..., where Br = 4cos2(7r/r) is the Tutte-Beraha number. The 
same caveat given above again applies; we have not made an exhaustive search for other 
even tinier sliver regions. The small regions R^ and R^ extend into the Re{q) < half-plane. 
In region Ri the dominant term is a root of the ninth-order equation with coefficient k^^^ 
which we denote X^^m, so that W = (Xa^riY^'^- The density of chromatic zeros is lowest on 
the portions of B that pass through g = 1 and q = and is intermediate on the portion 
that passes through g = 3. 



9 General Features 

Our exact results exhibit a number of features: 

• For a given width Ly, the Xiyj with coefficient ciyj = /t*-'^'', where 1 < d < Ly + 1, 
of which there are np{Ly,d), are roots of an algebraic equation of degree at most 
np(Ly,d). For example, for Ly = 3, the X^j with coefficients k^'^\ d = 4,3,2, and 1 
are, respectively, roots of equations of degree 1,3,5, and 4. This case illustrates the 
possibility that the degree of the corresponding equation is equal to np{Ly, d). A case 
illustrating the possibility that it is less is, e.g., that for Ly = 4, where the X^j for 
2 < J < 5 with coefficients n^^'^ are roots of a cubic equation and a linear equation. 

• Our general structural analysis shows that for a given Ly, there is a unique term Xl^^i 
with coefficient of maximal degree, and this coefficient is From our calculations 
we infer the general result that XLy,i = (— 1)"^". 

• Another inference concerns the set of terms Xr with coefficients There are 
Ly of these terms as shown in ( |3.2.8| ). Let X^^j = {—l)'^y~^^XLyj. Then the Xl^/s 



and hence the Xr ,'s with coefficients k^^^^ can be calculated as follows. Denote the 
equation whose solution is X^^j as f{Ly,^). Thus, /(1,0 = ^ ^ (O' ~ 2) and 

/(2,o = /(i,0(e-g + 3)-i. (9.1) 
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The Alj^j 's for higher values of Ly are then given by 

/(L„0 = /(i^,-l,0(e-g + 3)-/(L,-2,0 forL,>3. (9.2) 

We find that (i) if Lj^ = 1 mod 3, one of the ^Ly/s with coefficient k^^^^ is Xiyj = 
(_l)£'«+i(g_2); (ii) if = 2 mod 5, two others are {-l)^y+'^ times the roots of /(2,^); 
(iii) if Lj, = 3 mod 7, three others are (— 1)-^!/+^ times the roots of /(3, and so forth. 

The term that is dominant in region Ri is the (real) root of an equation with degree 
np{Ly, 1) = Ml^ of maximal magnitude and has coefficient 

The locus B for the ^ oo limit of the family GoiLy x L^) crosses the real axis 
and divides the q plane into various regions. This is in accordance with our earlier 
conjecture that a sufficient (not necessary) condition for families of graphs to yield 
a B that divides the q plane into different regions is that it contain global circuits, 
equivalent to periodic longitudinal boundary conditions for strip graphs. (Recall that 
for strip graphs with free longitudinal boundary conditions, B does not, in general, 



cross the real axis or divide the q plane into regions containing this axis p8|.) 



• In each case where we have an exact solution, we find that B crosses the real axis at 
q = I and q = 3. 

The finding that the locus B crosses the real axis on the left at g = 1 is different than, but 
related to, our earlier finding 0, 43], |45|-|49|, |5^, 53] that for the ^ oo limits of lattice 



strip graphs with periodic longitudinal boundary conditions and either free or transverse 
transverse boundary conditions, B crossed the real axis on the left at g = 0. For brevity of 
notation, we shall denote periodic and free longitudinal boundary conditions as PBC^. and 
FBCa;, respectively, and similarly, periodic and free transverse boundary conditions as PBCj, 
and FBCy. The present behavior can be understood as a consequence of the property that 
all of the vertices on the upper side of the strip are connected to a single external vertex. 
This is especially evident for the lowest case Ly = 1; here, as discussed in 0], it simply 
shifts the B for the ^ oo limit of the circuit graph, ]q' — 1] = 1, one unit to the right in the 
q plane to form the locus (|5.6|) . Hence, in particular, the crossing that would be present at 
g = for the circuit graph is shifted to a crossing at g = 1. This suggests that the crossings 
that are present at g = 1 for all of the cases for which we have exact solutions are universal 
in the same sense that the crossing on the left at g = was universal for the strip graphs 
with PBC^ and either FBCj, or PBC^. 

One of the most interesting features of the present work is the property that for all of 
the cases that we have studied, the point at which B crosses the real axis on the right, gc. 
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is also the same for different strip widths Ly. This is different from our ear her findings with 
lattice strip graphs with PBC^ and either FBC^, or PBCj^. There, we found that the generic 
behavior was, for a given set of these boundary conditions, and a given type of lattice, that 
qc{Ly) depended on Ly. For example, for strips of the square lattice with PBC^, and FBC^;, 
Qc = 2 for Ly = 1 and for = 2 [g, while ~ 2.34 for Ly = 3 ^ and Qc ~ 2.49 for 
Ly = 4 ||5^. From these results, we conjectured that Qc is a monotonically non-decreasing 
function of Ly for lattice strips with PBCx and FBC^ [|^. In contrast, we showed that Qc 
was not a monotonic function of Ly for lattice strips with PBC^; and PBC^. For example, 
gc = 2 for Ly = 2 (here the PBC^ and FBC^ results coincide), and qc then increases to the 



value Qc = 3 for Ly = 3, but decreases to qc — 2.78 for the next greater width, L„ = 4 [p3 
We inferred from these results that the crucial property (aside from having PBC^; so that 
B is guaranteed to cross the real axis and define a qc) is the type of transverse boundary 
conditions. The transverse boundary conditions for the Gd family of strip graphs are of 
neither pure free nor periodic type; the lower side of the strip has free boundary conditions 
while the vertices on the upper side are all connected to a single external vertex. We can 
characterize the strips with PBC^; and the various transverse boundary conditions by the 
number of sides with free boundary conditions; this number is 2 for FBCy, 1 for our present 
Go family, and for PBC^. It is therefore plausible that the above-mentioned conjecture 
that qc is a non-decreasing function of Ly for (PBCa;,FBCy) strips also holds for Gn strips. 
Accepting the validity of this conjecture, one can immediately infer the universality of the 
property qc = 3, independent of Ly, as follows. This property qc = 3 holds for at least 
one value of Ly, say, Ly = 1. One also knows that in the limit Ly ^ oo, qc = 3 for the 



(infinite) square lattice [£T1- Therefore, since the conjecture states that qc cannot decrease 
as Ly increases and since exact results show that it has already achieved its Ly = oo value 
at = 1, it must remain at this value for all larger widths Ly. 

As Ly increases toward infinity, one would expect that, in a certain sense, the effect of 
the single external vertex connected to each of the vertices on the upper side of the strip 
would become negligible. It would follow from this reasoning that in the limit Ly — > oo, the 
left-hand complex-conjugate arms of B would curve around, approaching the point q = 0, 
and finally meet at this point to form a closed component on B. Our results are consistent 
with this conjecture; we find that for the two largest values of Ly, where these left-hand 
arms are present (each with its own small region), they do approach the origin more closely 
for Ly = 4 than for Ly = 3. 
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10 Comparisons of W{Ly, q) for different L 



y 



As we have before for various infinite- length, finite lattice strips ^ , it is of interest to 

compare our exact solutions for IV at a given q, for various widths Ly, with the corresponding 
value of W for the infinite lattice, equivalent here to the limit Ly ^ oo. The range of q 
values for which this comparison is useful is determined by the values for which there is 
nonzero ground state entropy, W > 1 associated with the property that the number of ways 
of coloring the vertices of the strip such that no two adjacent vertices have the same color 



increases exponentially fast as as the number of vertices n ^ oo (cf. eq. |1.2|) . Here the 
limit n ^ oo is taken by sending oo. Although the chromatic number x alternates 

between 3 and 4 depending on whether is even or odd, one can take the limit of infinite 
length with even, and hence this sequence of families of graphs has x = 3^ the above- 
mentioned proper coloring can be carried out for g > 3. We thus make the comparison of the 
values of W{q) for various widths with the values of W{sq, q) for the infinite square lattice 
for this range q > 3. Except for W{sq,q = 3) = (4/3)^/2 the values W{sq,q) are not 
known exactly; we use the Monte Carlo calculations that were performed as part of |^3[| for 
this purpose. For the comparison, as before, we define the ratios 

^ H--(G.(L„xoo.,)) ^ 

W[sq,q) 

The results are shown in Table ^. 
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Tabic 5: Comparison of values of W{GD{Ly x oo),^) with W{sq,q) for 3 < g < 10. For each value of q, 
the quantities in the upper line are identified at the top and the quantities in the lower line are the values 
of Rw{Ly,q). 





W{l,q) 


W{2,q) 


W{3,q) 


W{4,q) 


W{sq,q) 


3 


1 


1.27202 


1.36138 


1.40596 


1.53960.. 




0.6495 


0.8262 


0.8842 


0.9132 


1 


4 


2 


2.16831 


2.22345 


2.25117 


2.3370(7) 




0.8558 


0.9278 


0.9514 


0.9632 


1 


5 


3 


3.12490 


3.16628 


3.18711 


3.2510(10) 




0.9228 


0.9612 


0.9739 


0.9803 


1 


6 


4 


4.09973 


4.13294 


4.14963 


4.2003(12) 




0.9523 


0.9761 


0.9840 


0.9879 


1 


7 


5 


5.08310 


5.11082 


5.12473 


5.1669(15) 




0.9677 


0.9838 


0.9891 


0.9918 


1 


8 


6 


6.07124 


6.09503 


6.10695 


6.1431(20) 




0.9767 


0.9883 


0.9922 


0.9941 


1 


9 


7 


7.06236 


7.08318 


7.09361 


7.1254(22) 




0.9824 


0.9911 


0.9941 


0.9955 


1 


10 


8 


8.05545 


8.07396 


8.08323 


8.1122(25) 




0.9862 


0.9930 


0.9953 


0.9964 


1 
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The property that all of the vertices on the upper surface of the strip, as represented in 
Fig. m are connected to an external vertex has the effect of restricting the coloring of the 
strip and thereby reducing the value of W^GoiLy x oc),q) relative to the value for a strip 
of the square lattice with the same periodic longitudinal boundary conditions, width, and q, 
but free transverse boundary conditions, W{Gs{Ly x oo), FBCy, PBCx, q). In a theorem 
was proved that for fixed g, W{Gs{Ly x oo), FBCy, PBCx, q) is a monotonically decreasing 
function of Ly which thus approaches the value of W{sq, q) for the infinite square lattice 
from above. Here, because of the effect of the coloring restriction due to the connections 
of all of the vertices on the upper surface to the single external vertex, the situation is 
different. Indeed, we find numerically that (for q > q^) over the range 1 < Ly < 4 where we 
have obtained exact solutions, W^GoiLy x oo),g) is an increasing function of Ly for fixed 
q in the range q > qc where . There are two effects depending on the strip width Ly that 
one can identify here, and these act in opposite directions. The first is that as the width 
increases, the restriction, per site, due to the connections of the vertices on the upper side to 
the external vertex is ameliorated, since these vertices occupy a smaller fraction of the total 
number of vertices as Ly increases. This is the dominant effect and tends to increase W. 
The second is that on the other (lower) side of the strip there is more freedom in coloring the 
vertices since they have no neighbors below them. These lower-side vertices occupy a smaller 
fraction of the total as Ly increases, and this tends to decrease W. It may be recalled that 
the proof given in for the monotonic decrease of W{Gs{Ly x oo),g) for fixed g > as 
a function of Ly for the case of free transverse boundary conditions relied upon this second 
effect. In that case, the vertices on both the upper and lower sides of the strip, which had 
greater freedom of coloring because of their reduced degree (number of neighboring vertices) 
occupied a commensurately smaller fraction of the total number of vertices as the strip 
width increased, and this yielded the result of the theorem. Here, however, one has a more 
complicated situation in which there are the two above-mentioned countervailing effects. 

As in previous works, we also record the value of W at certain other points in region Ri, 
in particular, g = and q = 1; these are given in Table Here we recall the order of limits 
used in our definition of W, ([17^ ). With the opposite order of limits, one has W{q)nq = 
for q = 0,1, and 2, for any Ly. 
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Table 6: Some additional values of values of \W{GD{Ly x oo),q)\. 



q 


W{l,q) 


W{2,q) 


W{3,q) 


W{A,q) 





2 


2.35829 


2.50361 


2.58283 


1 


1 


1.61803 


1.85334 


1.97672 



11 Conclusions 

Summarizing, in this paper we have presented a number of exact results for the partition 
function Z{Go, q,v) of the g-state Potts model and the T = antiferromagnetic special 
case given by the chromatic polynomial P{GD,q), on families of self-dual strip graphs Go 
of the square lattice of width Ly and arbitrarily great length with periodic longitudinal 
boundary conditions. We determined (i) the general coefficients cz,Ly,j and cp^Ly,d in terms of 
Chebyshev polynomials, (ii) the numbers nz{Ly, d) and np{Ly, d) of terms in Z and P with 
each type of coefficient, and (iii) the total number of terms Nz,Ly,\ and Np^Ly,x- Wc have 
pointed out interesting connections between the nz{Ly,d,) and Tcmperley-Lieb algebras and 
between the Np^Ly,x-, F = Z,P, and enumerations of directed lattice animals. We proceeded 
to present exact calculations of P for 2 < < 4 and to study the analytic structure of 
the resultant W functions in the complex q plane. A particularly interesting finding is a 
universal value, Qc — 3, for all of the widths Ly considered, which, furthermore, is equal to 
the value for the infinite square lattice. This is the first family of finite-width, infinite-length 
lattice strips for which we have found this type of universality of Qc- 
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12 Appendix 

In the text, we gave the A4J or the equations defining them for 1 < j < 5 for the Ly — A 
strip. In this appendix we give the equations for the terms X^j for 6 < j < 35. The A4J, 
6 < j < 14, are roots of a 9-degree equation 

- (g^ - llq^ + 51q^ - lUq + 101)^^ 
-(3g^ - 51q^ + 392q^ - 1756g^ + 4925?^ - 8597g^ + 8586g - 3750)^^ 
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-(3g^° - llq^ + 899g^ - 6307g^ + 29495g^ - 96135g^ + 221090g^ - 353821g^ 
+376297g2 - 23941 Ig + 68889)^^ 

- 40g^2 + 701g^^ - 7302g^° + 51090g^ - 255451g^ + 943367g^ - 2612569g^ 
+5434452g^ - 8391440g^ + 9348370?^ - 7109828g2 + 3304026g - 707059)^^ 
+ (4g^^ - 149g^^ + 2590g^3 - 27903g^2 + 208559g^^ - 1146435g^° + 4790250g^ 
-15497328?^ + 39142459?^ - 77176435g^ + 117779005?^ - 136551251?^ 
+116336276?^ - 68690109?^ + 25102423g - 4273553)^^ 
-(g - 2){<6q^^ - 226?^^ + 4004^^^ - 44291?^^ + 342389g^^ - 1961444g^^ 
+8613782?^° - 29579506?^^ + 80262537?^ - 172631424?^ + 293248190g^ 
-389119576?^ + 395172545?^ - 296710120?=^ + 155196867?^ - 50470243g + 7677301)^^ 
+ (g - 2)(4g^s - 167g^^ + 3293g^^ - 40754^^^ + 354781g^^ - 2307336?^^ 
+ 11619216g^2 - 46328775g^^ + 148278798g^° - 383788784?^ + 805173797g^ 
-1365777711g^ + 1858809772g'' - 2001697759g^ + 1667438591g^ 
-1036284383g^ + 452215858g^ - 123566922g + 15898678)^^ 
-{q - 2f{q^^ - 46g^^ + 994^^^ - 13427g^^ + 127231g^^ - 899290?^^ 
+4920180?^^ - 21334829?^^ + 74438607?^^ - 210900517?^° + 487345238?^ 
-918665623g^ + 1407005787?^ - 1735609937g^ + 1699139075?^ 
-1289913357?^ + 732092287?^ - 292275110?^ + 73185494? - 8642781)^ 
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-{q - 2y{q^^ - 34g^^ + 536g^^ - 5198?^^ + 34686?^^ - 168742^^° 

+618410g^ - 1738952g* + 3783595g^ - 6370774g^ + 8234267g^ - 8022902g^ 

+5703476g2 - 2792146g2 + 841318g - 117542) = . ( 
A4j, 15 < j < 26, are roots of a 12-degree equation 

+ {6q^ - 105g^ + 771g^ - 3057g^ + 6914^^ - 8454g + 4350)^^° 

+ (4g^ - 117g^ + 1470g^ - 10589g^ + 48615g^ - 148230g^ + 300934^3 

-392705g^ + 298868g - 100913)^^ 

+ {q^^ - 51q^^ + 1034g^° - 11831g^ + 87645g^ - 449620g^ + 1652070?^ 
-4404021g^ + 8481267g^ - 11530332?^ + 10515134g^ - 5776925g + 1445103)^^ 
-(6?^^ - 252?^^ + 4813^^2 - 55790g^^ + 440374?^° - 2511072?^^ + 10686844?^ 
-34528528?^ + 85175369g^ - 159716148g^ + 224144839?^ - 228278980?^ 
+159441815?^ - 68321288^+ 13538233)^^ 

+ (15g^^ - 637g^^ + 12667g^^ - 156625?^^ + 1348190?^^ - 8567954g^^ 
+41590251g^° - 157306151g^ + 468514489g* - 1102386000g^ + 2042068500g^ 
-2946106890g^ + 3244203521g^ - 2634966875g^ + 1487907358^^ 
-521546134g + 85418708)^^ 

-(20g^*^ - 910g^^ + 19575g^^ - 264460g^^ + 2514176g^^ - 17858194?^^ 

+98204545^^2 - 427439488?^^ + 1492563991g^° - 4211869748?^^ + 9626229067? 

30 



-17773020220?^ + 26306043996?^ - 30781790215?^ + 27840711649?^ 
-18773330488?^ + 8883458702?^ - 2630912469? + 366761469)^^ 
+ (15?20 - 749?^^ + 17749?^^ - 265398?^^ + 2808406?^^ - 22355292?^^ 
+ 138893776?^^ - 689698353?^^ + 2779911180?^^ - 9184226013?^^ 
+25005667770?^° - 56201345497?^ + 104079239038?^ - 157931128836?^ 
+194413953469?^ - 191128460116?^ + 146506854315?^ - 84366457483?^ 
+34322906908?^ - 8792343237? + 1066038269)^^ 

-(? - 2) (6?^^ - 324?^° + 8279?^^ - 133226?^^ + 1515528?^^ - 12965452?^* 
+86626388?^^ - 463261467?^^ + 2015823729?^=^ - 7215111460?^^ 
+21384133741?^^ - 52649552494?^° + 107694098636?^ - 182459513387?^ 
+254392095608?^ - 288780678892?^ + 262603966158?^ - 186683744725?^ 
+99924573125?^ - 37861221344?^ + 9048442752? - 1025103882)^2 
+ (? - 2)(?23 - 65?^^ + 1965?^^ - 36996?^° + 488991?^^ - 4842151?^^ 
+37397267?^^ - 231374113?^^ + 1168070271?^^ - 4873964424?^^ 
+16957771126?^^ - 49471796934?^^ + 121369147461?^^ - 250500702268?^ 
+434093944035?^ - 628647524990?^ + 754930635327?^ - 743026096643?^ 
+589226915539?^ - 367162228065?^ + 173056035863?^ - 57980197360?^ 
+12299658707? - 1241302954)^^ 
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+{q - 2f{q - 3)(2g22 _ iQg^2i ^ 2662?^° - 42156?^'^ + 472587g^^ 
-3990653g^^ + 26367347?^^ - 139766257g^^ + 604533952g^^ 
-2158390643g^^ + 6409094853^^^ _ I5896051134g^^ + 32980281410g^° 
-57170276403g^ + 82469188225g^ - 98281609520g^ + 95679352146g^ 
-74828758316^^ + 45863717844g^ - 21212341252^3 + 6958799896?^ 
-1442606958? + 142022490)^ 

+ (g - 2)^(^2 - 5? + 5)(g^° - 46g^^ + 999?^^ - 13621g^^ + 130785g^^ 

-940159?^^ + 5250952?^^ - 23336758?^^ + 83832989g^^ - 245863743?^^ 

+591990845g^° - 1172460575?'^ + 1906935986?^ - 2533386002?^ 

+2722446125?^ - 2330163861?^ + 1551214058?^ - 774015177?^ 

+272288617?^ - 60200487? + 6289073) = . (12.2) 
Finally, the A4J, 27 < j < 35, are roots of another 9-degree equation 

- (6?2 - 33? + 48)^^ + (15?^ - 168?=^ + 707?^ - 1327? + 933)^^ 

-(20?^ - 345?^ + 2450?^ - 9190?^ + 19218?^ - 21235? + 9659)^^ 

+ (15?^ - 360?^ + 3693?^ - 21232?^ + 75012?^ - 167016?^ + 229025?^ 

-176806? + 58742)^^ 

-(6?^° - 195?^ + 2732?^ - 21931?^ + 112366?^ - 385466?^ + 899022?^ 
-1410181?^ + 1425238?^ - 838361? + 217832)^^ 
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+{q - 3)(g^^ - 45g^° + 806?^ - 7979g^ + 49672g^ - 206895?^ + 59320 V 
-1177217g^ + 1590681g2 - 1397396?^ + 719432g - 164543)^^ 
+{q - 2){3q^^ - 108g^^ + 1726g^° - 16260g^ + 100858g^ - 434873g^ + 1338650g^ 
-2967964g^ + 4708729g^ - 5217472g3 + 3834660g^ - 1678782g + 330985)^^ 
+ {q -2){q- 3)(3g^2 - 93g" + 1295g^° - 10716g^ + 58721g^ - 224620g^ 
+615435g^ - 1217860?^ + 1728820^^ - 1718182?^ + 1135486?^ - 448175g + 79896)^ 
+(q - 2f{q^^ - 32q^^ + 465g^^ - 4062?^° + 23798?^ - 98757g^ 
+298728?^ - 667187g^ + 1100852?^ - 1326203?^ + 1134621?^ 

-653205?^ + 226903? - 35923) = . (12.3) 
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